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We show that the stochastic dynamics of a large class of one-dimensional inter- 
acting particle systems may be presented by integrable quantum spin Hamiltoni- 
ans. Generalizing earlier work jL], ^] we present an alternative description of these 
processes in terms of a time-dependent operator algebra with quadratic relations. 
These relations generate the Bethe ansatz equations for the spectrum and turn the 
calculation of time-dependent expectation values into the problem of either finding 
representations of this algebra or of solving functional equations for the initial val- 
ues of the operators. We use both strategies for the study of two specific models: 
(i) We construct a two-dimensional time-dependent representation of the algebra 
for the symmetric exclusion process with open boundary conditions. In this way 
we obtain new results on the dynamics of this system and on the eigenvectors and 
eigenvalues of the corresponding quantum spin chain, which is the isotropic Heisen- 
berg ferromagnet with non-diagonal, symmetry-breaking boundary fields, (ii) We 
consider the non-equilibrium spin relaxation of Ising spins with zero-temperature 
Glauber dynamics and an additional coupling to an infinite-temperature heat bath 
with Kawasaki dynamics. We solve the functional equations arising from the al- 
gebraic description and show non-perturbatively on the level of all finite-order 
correlation functions that the coupling to the infinite-temperature heat bath does 
not change the late-time behaviour of the zero-temperature process. The associ- 
ated quantum chain is a non-hermitian anisotropic Heisenberg chain related to the 
seven-vertex model. 
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1 Stochastic Dynamics and Quantum Systems 



One-dimensional stochastic reaction-diffusion processes are of both theoret- 
ical and experimental interest in a very wide context. They are well-known 
models both for reaction-diffusion mechanisms in physics and chemistry and 
for stochastic spin flip dynamics ||. More recently they have appeared 
through various mappings also as models for traffic flow ||, [6|, [JJ |8[], the 



kinetics of biopolymerization M, 10], reptation of DNA in gels O, 13 



14] , [T5| , PI , interface growth [[IT], |TS|| , diffusion in zeolites [jT^, ^(J and many 



other phenomena. Even in relatively simple models of driven diffusion such 
as the asymmetric exclusion process with open boundaries one finds a very 
rich dynamical behaviour involving dynamical and non-equilibrium phase 
transitions of various kinds [21, 22, 23, 24]. Exact solutions [§, ^| allow for 



a detailed understanding of cooperative phenomena in these classical many- 
body systems and provide insight in the role of inefficient diffusive mixing in 
diffusion-limited chemical reactions, in the dynamics of shocks and in other 
fundamental mechanisms which determine the behaviour of low- dimensional 
systems far from thermal equilibrium. 

A convenient and much used description of stochastic processes is in terms 
of a master equation for the probability distribution f(n; t) of the stochastic 
variables n. These variables represent the states in which the system may 
be found at any given instant of time. The master equation encodes the 
transition probabilities p(n' — > n) of moving from one state nf to another state 
n in one time step: the probability distribution f(n;t + At) = J2 n ' ~* 
R)f{BL',t) is just the sum of probabilities of finding the system in state nf_ 
at time t times the respective transition probabilities — > n). Thus the 
master equation expresses the probability of finding the system at time t + At 
in a given configuration n in terms of the probability distribution at time 
t. Such processes are Markov processes which may be constructed for the 
description of interacting particle systems |26[ . 



Since the master equation is linear in the probability distribution, it can 
be expressed as a vector equation in a "quantum Hamiltonian formalism" 
by mapping each state n of the system to a basis vector | n ) in a suitable 
vector space X. Thus the probability distribution becomes a vector | f(t) ) = 
Y^n f{n',t)\n) an d the master equation for a continuous-time process (At — > 
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0) takes the form of a Schrodinger equation in imaginary time 



j t \f(t)) = -H\f(t)). (1) 

For various interesting many-body systems the generator H of the (classical) 
stochastic time evolution turns out to be identical to the quantum Hamilto- 
nian of spin chains known from condensed matter physics |^7|, |28|, ^9], |3(], [31], 
|32|1 (see next section). The ground state of a stochastic Hamiltonian (which 
by construction has energy 0) corresponds to the stationary probability dis- 
tribution of the stochastic dynamics which is reached asymptotically at the 
late stages of the time evolution. Through this somewhat surprising connec- 
tion to quantum spin systems the stochastic dynamics become tractable with 
the tool box of quantum mechanics and progress may be achieved.^] 

A paradigmatic example of this relationship is the symmetric exclusion 
process [|6|, [37f1 . In this lattice model particles hop between lattice sites k, I 
with rates p^i = Pi t k- They interact via a hard-core repulsion which prevents 
the occupation of a lattice site by more than one particle. The stochastic 
quantum Hamiltonian obtained for this system |28| is the Hamiltonian for 
the isotropic spin- (1/2) Heisenberg ferromagnet 

H = -\T, Vki [°iof + a y k af + otf - 1] . (2) 
z k,i 

The manifest S77(2)-symmetry of H (which is not recognizable in the original 
master equation) can be used for obtaining exact results ^9], [38|. Moreover, 
in one dimension, the system with nearest neighbour hopping is integrable 
and can be solved by the Bethe ansatz P3I . 



The integrability is not a special feature of the symmetric exclusion pro- 
cess alone. Driven diffusion in one dimension {pk,k+i/Pk+i,k — Q 1) is 
described by the XXZ quantum spin chain which differs from (0) by an 



1 For discrete time dynamics the vector form of the master equation reads | f(t + At) ) = 
T\ f{£) ). For many interesting systems in one dimension T is the transfer matrix of a two- 
dimensional vertex model ||| |35|, [56| . The frequently- used notion "quantum Hamilto- 
nian formalism" is somewhat misleading in so far as for many stochastic applications the 
transition rates result in coupling constants which make H non-hermitian. Further, since 
/ is a (real- valued) probability and not a (complex) probability amplitude, the expectation 
values of the stochastic process are not the expectation values normally calculated in a 
quantum mechanical problem. 
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anisotropy A = (q + g~ 1 )/2 in the z — z coupling of the spins [p9| . Us- 
ing the Bethe ansatz and related methods many new exact results par- 
ticularly on the dynamical properties of the process have been obtained 
|1|, ||, m ||, ||, |6j. Moreover, it turns out that a 10-parameter 
class of reaction-diffusion systems of identical particles |J2| and various sys- 



tems of non-identical particles ^] are described by integrable quantum 



chains. Unfortunately however, the Bethe ansatz requires knowledge of some 
reference eigenstate of H and integrability has so far failed to provide con- 
structive methods to calculate even the ground state of spin chains with 
boundary fields that destroy the reference state. Such boundary fields are of 
importance in stochastic dynamics of non-equilibrium systems as they allow 
for a modelling of open systems which are in contact with external parti- 
cle reservoirs at their boundaries and thus allow for maintaining a particle 
current through the system. 

In a very different approach the ground states of one-dimensional spin 



Hamiltonians are formulated in terms of matrix product states [48], f|9|, |K| pi 
where the ground state wave function is expressed in terms of a trace over 
a product of matrices. They may be seen as representations of an operator 
algebra which is determined by the requirement that by acting with the 
Hamiltonian on this state one obtains an eigenstate (the ground state) of H . 
Applied to stochastic Hamiltonians one obtains in this way the stationary 
distribution of a stochastic process |52|, |53|, |54|, [55], [56], |57jl.p] 



By constructing an infinite-dimensional representation of the stationary 
matrix algebra for the asymmetric exclusion process with open boundaries 
Derrida et al. J52] produced the same exact results that were obtained in- 



dependently by Schiitz and Domany |E| using a different method. In fact 



with hindsight our treatment may be seen as a representation-free solution of 
recursion relations that one can derive from the matrix algebra. It is the aim 
of this work to apply both strategies, viz. (i) construction of a matrix repre- 
sentation or (ii) solution of equations resulting from the algebraic relations 
alone, to a matrix product treatment of the dynamics of reaction-diffusion 
systems. Specifically we consider the symmetric exclusion process with open 
boundaries and a reaction-diffusion mechanism which is equivalent to a spin 
relaxation model. 



2 For non-periodic systems with boundary fields one does not take a trace, but a suitably 
chosen scalar product [p2| . 
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The extension of the stationary matrix approach of Derrida et al. to a dy- 
namical description requires one new idea. This is the introduction of auxil- 
iary matrices S, T which do not appear in the calculation of expectation 
values, but are necessary to formulate the dynamical algebra which is deter- 
mined such that the probability distribution satisfies the full time-dependent 
master equation. The special case of this construction where S + T = 
allowed among other things for the rederivation || of the Bethe ansatz equa- 
tions |58[ for the spectrum in the symmetric case with open boundaries and 
for the rederivation []59 of the spectrum in the asymmetric case with periodic 
boundary conditions [|3|. With the auxiliary matrices as additional ingre- 



dient the extension of the dynamical matrix product approach to reaction- 
diffusion systems becomes straightforward |5U| (see below). 

I will first show (Sec. 2) how reaction-diffusion systems of identical hard- 
core particles are related to a generalized Heisenberg chain. Its spectrum 
can be obtained from the Bethe ansatz. This is a simplified rederivation 
of some results obtained earlier fl32"| . Then in (Sec. 3) I will generalize the 



operator approach to the general reaction-diffusion problem of identical hard- 
core particles with nearest neighbour interaction in one dimension. As an 
application I will return to the symmetric exclusion process and present a 
two-dimensional representation of the time- dependent operator algebra.0 In 
Sec. 4 I will solve a non-equilibrium spin relaxation model introduced by Droz 
et al. |)T| by solving functional equations arising from the algebraic relations 
of the time-dependent matrix algebra. This treatment does not require the 
calculation of representations of the algebra. In Sec. 5 the main results are 
summarized and some open questions are pointed out. 



2 Integrable Reaction-Diffusion Processes 

We will consider stochastic reaction-diffusion processes of identical particles 
with hard-core repulsion moving on a ring with L sites. Even though part 
of our approach generalizes to arbitrary lattices [3~2~] we will study here only 
one-dimensional systems with nearest neighbour interaction. The stochastic 
variables of the system are the occupation numbers n = {n^} where nk — 0,1 



3 Sections 2 and 3 of this paper are not really new. They constitute the bulk of the paper 
]60| ] which was presented at the Satellite Meeting to Statphys 19 on Statistical Models, 
Yang-Baxter Equation and Related Topics, at Nankai University, Tianjin (August 1995). 
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Table 1: Bulk reaction and diffusion rates for nearest neighbour exclusion 
processes of identical particles. The numbers the rate of change of 

the occupation numbers {n k ,rik+i}- 



Process 


Rate 




Process 


Rate 




01 -> 10 




diffusion 


10 -> 01 


<*23 


diffusion 


11 -> 01 


024 


coagulation 


01 -> 11 


042 


decoagulation 


11 -> 10 


034 


coagulation 


10 -> 11 


043 


decoagulation 


00 -> 01 


a 2 i 


creation 


01 -> 00 


Ol2 


annihilation 


00 -> 10 


a 3 i 


creation 


10 -> 00 




annihilation 


00 -> 11 


a 4 i 


pair creation 


11 -> 00 


Ol4 


pair annihilation 



indicates whether site 1 < k < L in the lattice is occupied or empty. At a 
given time t the state of the system is completely described by the probability 
distribution /(n; t). In this class of models there are ten possible reactions in 
addition to right and left hopping (diffusion), so altogether one has to specify 
12 independent rates > (Tab. 1). 

The stochastic dynamics are defined by the master equation 

4/(ffif) = E Hmn!)f(nf;t) - w(nf;n)f(n;t)\ (3) 
dt ~, 

where the reaction-diffusion rates w(n;n/) for a change from configuration 
n/ — > n are equal to the sum 
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fc=l 



a 12^n k ,0^n k+1 ,0 + a 32^n k ,l^n k+1 ,0 + a 42<^n fc ,l^n fc+1 ,l 
a 13Vn k fifin k+1 ,0 + a 23^n k ,0^n k+1 ,l + a 43^n fc ,l^n fc+1 ,l 
(lu8n h ,o8n h+1 ,0 + a 240n k flOn k+1 ,l + O34^n fc ,l5n fc+1 ,0 



5 n'„0^ +1 ,l 

<k,i<k +1 ,o 



+ 

+ 
+ 

}■ 



This somewhat lengthy expression becomes more compact in the quantum 
Hamiltonian formalism ([I]): To each configuration n a vector | n ) which, to- 
gether with the transposed vectors (n |, form an orthonormal basis of (C 2 )® L . 
In spin language this corresponds to a mapping to a spin 1/2 chain by identi- 
fying a vacancy (particle) at site k with spin up (down) at this site. The prob- 
ability distribution is then given by the state vector | f(t) ) = J2n /fe 1)\ n) 
and the formal solution of the master equation (Q) in terms of the initial 
distribution | /(0) ) is given by | f{t) ) = exp(—Ht)\f(0)). The stochastic 
dynamics are defined by the master equation (|l|) with|3l 



H 



L 
k=l 



(4) 



where the matrices hk act non-trivially only on sites k, k + 1 and are given 

by 



hk — — 



( an 

021 
031 

V a 4 i 



Ol2 
«22 
«32 

a 42 



«13 
«23 
«33 

a 43 



a 14 ^ 

«24 
«34 

a 44 / 



(5) 



with a 



The connection of H to the Heisenberg quantum chain becomes apparent 
by the similarity transformation H = QVHV" 1 ^ 1 with V = exp(5' + ) where 
S + = J2k=i s t an d s fc = ( a k i zcr fe)/2 are the spin lowering and raising 
operators acting on site k and with $ = exp(£ J2k ^ a l) where £ is a suitably 
chosen constant |32|. On the ten parameter submanifold defined by 



a 34 — (221 + a 41 + a 12 + a 32 — a 23 
«24 = Ct31 + ^41 + <3l3 + a 23 — &32 

the transition matrices have now the structure hk 



a 43 — a lA 
O42 — a 14- 



(6) 
(7) 

+ /i^T. Here h? xz 



commutes with S z = Y.k=i and h k is a sum of two parts which lower the 
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^-component of the spin on sites k, k + 1 by one and two units respectively. 
So one finds 

H = H xxz + H (8) 

where H xxz is the Hamiltonian of the anisotropic Heisenberg ferromagnet 
with twisted boundary conditions in a magnetic field. The crucial observation 
is that H~ does not change the spectrum ofH xxz , since H xxz may be block- 
diagonalized into blocks with fixed quantum number S z and H connects 
only blocks of given S z with blocks with quantum numbers S z — 1 and S z — 2.0 
Quantities of interest are expectation values (i.e. r-point correlation func- 
tions) ( n kl (t) . . . n kr (t) ) fo = (s \n kl . . . n kr e~ Ht \ /(0) ) which give the proba- 
bility of finding particles on the set of sites {k\, . . . , k r } at time t, if the initial 
distribution at time t — was f . Here ( s \ — X) n (n| and n k = (1 — of)/ 2 
is the projector on states with a particle on site k. ^From the Bethe ansatz 
one finds now that the spectrum has an energy gap (i.e. inverse correlation 
time) /j/ = 4a 41 + 2(a 2 i + 031) + 012 + ^13 — «42 — a 43 > 0. If // = the 
dynamical exponent turns out to be z = 2. Note also that V transforms a 
r-point density correlation function into a matrix element in the sector with 
r down spins. Since H~ only creates down spins, only transformed initial 
states with I < r down spins will contribute to the correlation function. This 
surprising simplification allows for an exact calculation of the local average 
density for any initial state even though we are dealing with a non-trivial 
interacting many particle system . 



3 The Dynamic Matrix Ansatz 

The results of the last section involve the constraints (|]), (0) and do not ap- 
ply e.g. for the asymmetric exclusion process. Also this model is integrable, 
but a calculation of time- dependent correlation functions has not yet been 
achieved. In order to solve this problem we now formulate a dynamic matrix 
ansatz for the general react ion- diffusion system defined by (|j) and (|5p, gen- 
eralizing earlier work Jl|, Q] for diffusion only. Instead of periodic boundary 
conditions we consider a system with open boundaries where particles are 
injected (absorbed) at site 1 with rate a (7) and at site L with rate S {(5). 

4 This mechanism was first noticed in a similar context in Alcaraz et al. |p2|. 
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bi 



I hk with suitably chosen boundary matrices 



Therefore H 

The ansatz is to take | f(t) ) = (( W \{Uk=i W) + )}| )l V))/Z L 

where | 0) is the state with all spins up and D,E are time-dependent ma- 
trices satisfying an algebra obtained from the master equation (JJ). The 
(time-independent) vectors ((W| and \V )) on which D and E act are 
determined from the boundary terms in the master equation and Zl = 
(( W \C L \ V )) where C = D + E is a normalization. In this framework 
the r-point density correlation function is given by ( if) . . . n kr (t) ) / = 
(( W \C kl - 1 DC ka - kl - 1 D . . . C L - kr \V))/Z L . Therefore, given a matrix repre- 
sentation of the algebra satisfied by D, E, the computation of time-dependent 
correlation functions is reduced to the much simpler calculation of matrix el- 
ements of a product of L matrices. 

It is easy to see that (|I|) is solved if for each pair of sites one satisfies 



1 d 
2~dl 

(S + Ta k )(E + Da k+1 



+ h k )(E + Da k )(E + Da k+1 )\0) = 

(E + Da k )(S + Ta k+1 )]\0) (9) 



where S, T are auxiliary operators satisfying 



((W\ 

.1 d 
2~dt 



(l± + bl )(E + Da^ + (S + Ta^\0) 



+ b L )(E + DaZ)-(S + Taj:)\0) 



V)) 




0. 



(10) 

(11) 



By comparing each of the four terms in @ proportional to | ) , a k \0), 
a fe+1 |0) and a k o~ k+1 \0) resp. one obtains four quadratic relations for the 
operators D,E,S,T. Eqs. ( |T0D and (|TT|) give two pairs of equations which 
define ((V^l and | V)). Introducing 



.4 
B 
B 
A 



(i) 
j 

(i) 
i 

(2) 

i 

(2) 



-(o 2 i + a 31 + a 41 )E 2 + a 12 ED + a 13 DE + a 14 D 2 (12) 

a 2 iE 2 - (012 + a 32 + a 42 )ED + a 23 DE + a 24 D 2 (13) 

a 31 E 2 + a 32 ED - (a 13 + a 23 + a 43 )DE + a 34 D 2 (14) 

a 4 \E 2 + a 42 ED + a 43 DE — (ai 4 + a 24 + a 34 )D 2 . (15) 
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one finds 



~E 2 -[S,E] = A« (16) 

-^-ED-SD + ET = B {1) (17) 

-^DE -TE + DS = 5 (2) (18) 

2dt K ' 

\j t D 2 -[T,D] = (19) 



and 



{{W\\^j t E-aE + 1 D + S^ = (20) 

(( W \[\j t D + aE ~^ D + T ] = ( 21 ) 
^E-6E + (lD-s}\V)) = (22) 

{~D + 6E-PD-t\\V)) = 0. (23) 

One may reduce this algebra by assuming that C is time-independent and 
has a representation where it is invertible. Eqs. © then imply [C, S+T] — 



and ([HI,© imply (( W |(5 + T) = = (S + T)| V )). This can be solved 
by assuming S + T = 0, which, as I would like to stress, is not the most 
general choice. Now one can express S in terms of C and D and is left 
with only two further relations to be satisfied by D and C and two relations 
defining (( W \ and | V)). In particular, if (0) and (§) are satisfied, there is 
one relation involving D which is linear in D and one relation quadratic in 
D. For the symmetric exclusion model 023 = 032 = 1/2 this dynamic al- 
gebra yields eigenvalue equations for the corresponding XXX-Hamiltonian 
with integrable, but non-diagonal, symmetry breaking boundary fields]!], [J. 
However, no matrix representation has been found yet. This raises the ques- 
tion whether non-trivial representations do exist at all. 

As I will show here for the first time, the answer to this question is yes, at 
least with some restrictions on the injection and absorption rates. Choosing 
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a basis where C is diagonal one finds the representation 





with e = (a + /? + 7 + S)/2, c=l-a-7=(l-/3- S)~ x , d = a/(a + 
7 ) = 5/(0 + 5) and ((W|, arbitrary but ({W\V)) 0. In this 

representation A is an arbitrary parameter specifying the initial distribution. 
One may also use it for the construction of (right) eigenstates of H, since the 
expression (( W \E kl ~ 1 DE k2 ~ kl ~ 1 . . . E L ~ kr \ V )) is a superposition of wave 
functions ty €i (ki, . . . ,k r ) of eigenstates with eigenvalues £j. The argument is 
the position of r down spins on sites k\, . . . , k r . Taking A = corresponds 
to taking the stationary distribution as initial state. This is an eigenstate 
with energy 0. The terms proportional to A give the wave function for an 
eigenstate with energy e. The quantity 1/ (In |c|) plays the role of a spatial 
correlation length. 



4 Non-equilibrium spin relaxation 

A phenomenon of wide interest in physics and chemistry is the growth of 
domains in non-equilibrium two-phase systems. The best-known example is 
perhaps the Ising model with domains of up- and down spins, separated by 
domain walls. The energy of the Ising model is given by the nearest neighbour 
sum E = —JJ2 s i s j- Since the creation of a local domain wall costs an energy 
J the system tries to organize itself at low temperature into large domains 
of uniform magnetization. Starting from a high-temperature equilibrium 
state with many domain walls and quenching to low temperatures leads to a 
coarsening process: Small domains of uniform magnetization merge to form 
larger domains since then the total length of the domain walls and thus the 
energy decreases. 

Glauber |63[ introduced spin-flip dynamics which ensure that the system 



reaches the equilibrium distribution at temperature T = 1/(3 of the one- 
dimensional Ising model. In this model a spin within a domain of equal 
magnetization is flipped with a rate fi = 1 — tanh/3J, whereas a spin in a 
region of opposite magnetization is flipped with a rate A = 1 + tanh/SJ. At 
domain boundaries spins are flipped with unit rate, since no change in energy 
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is involved. This process can be visualized in the following way: 



T T T —> T I T and I I I I t i w ^h rate /i 
t I T -> T T T and | t I -»• I I I with rate A 
T T I - T I I and I | t - | T T with rate 1 

Glauber dynamics can also be seen as a reaction-diffusion system. One 
simply identifies an up-spin with a vacancy and a down-spin with a particle. 
In one dimension at zero temperature the process can then be described as 
follows: 

A or A -> A A with rate 1 
A or A -> 00 with rate 1 



This can obtained by a translational rearrangement of the three-site interac- 
tions in terms of two-site processes. One realizes then that Glauber dynam- 
ics can be represented by a stochastic Hamiltonian of the form ([5]) with 
«i2 = a i3 — a 42 = 043 = 1- Furthermore, these rates satisfy the constraints 
(|7p. We stress that this relation to a reaction-diffusion system is not really 
a mapping, but just a certain choice of language which we use in order to 
make contact with Sections 2 and 3. There are two different non-trivial map- 



pings 1 64, 55] to the process of diffusion-limited annihilation which has little 



in common with the process described here. These mappings are useful as 
they show that Glauber dynamics can be described and solved in terms of 



free fermions [27, QQ, |>7], p{J and that the associated quantum chain 



is a non-hermitian anisotropic Heisenberg chain related to the seven-vertex 
model. 

The Glauber relaxation rules involve single spin flips and thus do not 
conserve the total magnetization. Kawasaki [nj introduced spin-exchange 
dynamics which also lead to an equilibrium Ising distribution, but which do 
conserve the total magnetization. At infinite temperature, these dynamics 
reduce to simple exchange of neighbouring spins with some rate a 2 3 = a 32 = 
(, i.e. to the symmetric exclusion process described above. Due to the lack 
general theorems on the dynamics of non-equilibrium systems it is now of 
interest to investigate a spin system the behaviour of which results from a 
coupling to two heat baths at different temperatures - one leading to zero- 
temperature Glauber dynamics, the other to infinite-temperature Kawasaki 
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dynamics [|6T]. In such a situation there is a competition: The diffusion 
process tries to disorder the system, while the Glauber process tries to create 
an ordered system of uniform structure. Hence the questions arise, which 
process wins, and how is the stationary state reached. 

In one dimension this problem was addressed by studying the dynami- 
cal spin-spin (= particle-particle) correlations for a translationally invariant 
initial state ||61|| , using the fact that the equations of motion for correlation 
functions decouple into closed subsets. By solving the equations Droz et al. 
could show that at any (finite) value ( of the coupling strength the system 
orders and that the spin-spin correlation function behaves at large times like 
the zero-temperature Glauber correlator. Here we use the matrix product 
ansatz to prove non-perturbatively that this remains true for all correlation 
functions of finite order. The leading contribution to time-dependent corre- 
lation functions (for large times) is always the zero-temperature Glauber cor- 
relation function |7T| . Corrections resulting from the coupling to the infinite- 
temperature heat bath are of subleading order 1 / \ft (relative to the leading 
contribution) .f\ 

To prove our assertion we consider the matrix algebra describing the 
process. Since we are more interested in spin variables we introduce D = 
C — 2D = E — D. We restrict ourselves to translationally invariant initial 
states. Because of translational invariance time-dependent spin expectation 
values are given by a trace over matrices 

(<W---<W)/o = Tr {C k ^ 1 DC k ^- 1 D...C L -^}/Z L (25) 

where = Tr C L . We reduce the algebra (16) - (p~9|) as in the case of the 
symmetric exclusion process by setting S+T = 0. Eliminating S leads then to 
the algebra generated by C, D and C^ 1 with the relations CC^ 1 = C~ 1 C = 1, 

5 We mention in passing that this process is a simple toy model of growing tissue cell 
populations [J72|. The decoagulation process A0,0A — > A A with unit rate describes cell 
division, while the particle hopping with rate C corresponds to the diffusive motion of cells 
in their environment. In addition to that we allow for a death process A0, %A — > 00 with 
rate q which kills both the original cell and its offspring during the decoagulation (cell 
mitosis) . It is intuitively clear that for q < 1 the cell population will grow until all space is 
covered, while for q > 1 the population will eventually die out. Therefore it is of interest 
to study the case q = 1 when creation of offsprings and the death process balance each 
other. The process leads to an ordered state also in three dimensions [ |F3| . This implies 
that cither all tumor cells die, or, with equal probability, cover the whole available space. 
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d/(dt)C = and 

jD = (l + (CDC- 1 + C~ x bC - 2D) (26) 
2(1 - A) = DC^DC- 1 + C- 1 DC- 1 D -2AC- 1 D 2 C- 1 (27) 

with the constant A = £/(l + £) determined by coupling ratio £. The r-point 
spin correlation function can now be rewritten 

( < W • • • *t («) ) = Tr {£ fcl . . . D^C^/Zr, (28) 

where Dk = C k ~ 1 DC~ k and fcj+i > fcj. These relations provide an alterna- 
tive, purely algebraic definition of the spin relaxation process. 

Relation ( p6|) is linear in D and we procede by constructing the Fourier 
transforms T> p = Y^k e tpk Dk to reformulate the algebra in terms of the Fourier 
components. Since 

CVpC- 1 = e~ ip V p , (29) 
the time-dependence of T> p is now simply obtained from (|26|) 



V p {t) = e -^P p (0) (30) 

in terms of the initial matrix T> p (0) and the "energy" 

e p = 2(1 + 0(1- cos p). (31) 

/From (g7p follows 2(1 - A)£(p) = J dp'V p ,V p _ p ,{l + e" ip - 2Ae 4p '-* p ). 
Since this relation holds for all times (and all p), the integral can be divided 
into separate time-components, each of which must vanish. If p\ and P2 are 
non-zero this leads to 

V P1 V P2 = S (p 2 , pi ) V P2 V Pl ( 32 ) 

with the two-body scattering matrix 

X _l_ e ipi+ip2 _ 2Ae ip2 



known from the usual anisotropic Heisenberg chain |4(J 



This relation for T> p is derived for the dynamic components with p ^ 0. 
Hence the static term in the l.h.s. of ( |27|) does not reappear in (0), but 
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in a different relation involving the static parts T>q = J2 n C n ~ 1 DC~ n and 
X = J2n n C n ~ 1 DC~ n . These quantities need separate treatment in a similar 
way as the static components of the operators for the symmetric exclusion 
process with open boundaries [[25], 74[ . However, since the stationary state 



(all spins up or all spins down) is not interesting we do not consider the static 
Fourier components. 

The momentum space formulation of the algebra provides another equiv- 
alent formulation of the process. To calculate expectation values we do not 
search for a representation of this algebra but use a different strategy. In 
terms of the Fourier components the correlator (|28|) reads 



<<(*)•••<(*)> = (ft/ f*e-**«-**) T({ft}) (34) 

where the so far undetermined matrix element 

T({ Pi }) = Tr {XV(0) . . .V P J0)C L }/Z L (35) 

depends only on the initial distribution. 

First consider the one-point function ( cr| (t) ) . Because of translational 
invariance, this local magnetization is independent of space. This is reflected 
in the invariance of the trace under cyclic permutation. Thus (^) for the one- 
point function together with ( |2"9"D imply e ip = e 2ip = e 3ip = . . . = e ipL = 1 
and therefore T(p) = 2nc5{p). This yields ( cr|(t) ) = c for all times and 
restates nothing but the known result that the average magnetization remains 
constant under the time evolution of the spin relaxation model pi]]. The 
initial magnetization m fixes c = m . 

We obtain a non-trivial result for the two-point correlator. We use ([29]), 



(33) and the cyclic property of the trace and find 



= S(p 2 , Pl )Tr{V P2 V Pl C L }/Z L 

= e^ L S(p 2 , Pl )Tr {V P2 C L V P1 }/Z L 

= e ip ^ L S(p 2 ,p 1 )T( Pl ,p 2 ). (36) 

A similar cycling procedure gives a second relation 

T(pi,p 2 ) = e ip2L S( Pl ,p 2 )T( Pl ,p 2 ). (37) 
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Using S(pi,p 2 ) = S 1 (p 2 ,pi) we conclude that T( P i,p 2 ) is non-zero only if 
the momenta p± )2 satisfy the Bethe-ansatz equations JT(| 



e ipiL = S( Pl ,p 2 ) (38) 
e ip2L = S(p 2 , Pl ). (39) 

Translational invariance requires also Pi + p 2 = 0. 

Moreover, the first of the equations (|36| ) yields a functional equation for 
the matrix element 

T(p 1 ,p 2 ) = S(p 2 ,p 1 )T(p 2 ,p 1 ). (40) 
This functional equation for T(pi,p 2 ) is solved by the Bethe wave functions 



y piiP2 (h,h) = A(h,l 2 ) ( e ^ 1+ ^ 2 + S(p 2 , Pl )e^ h+ip ^) , (41) 

for k 2 > k\ with some amplitude A(li,l 2 ) determined by the initial value of 
the correlator. It can be calculated from the integral representation ( |34]) of 
the full time-dependent correlator by setting t = 0. 

For a finite system the integral has to be replaced by a sum over the 
solutions of the Bethe ansatz equations (|38|), (|39|) . However, in an infinite 
system the set of solutions becomes dense. The only subtlety arises then 
from the bound states defined by the pole of the scattering amplitude S. 
This pole corresponds to the two-particle bound state already known from 
the original solution of Bethe J^J for A = 1. One can fix the contour of 
integration by setting A(li,l 2 ) = (71^(0)77,^(0) ). This gives 

T( Pl , P2 ) = £ (n h (0)niM ) (e ip ^ +ip2h + S(p 2 , Pl )e^ h+ip ^) (42) 

with the sum being restricted to the domain l 2 > li. We prescribe the 
appropriate contour of integration by isolating in S the constant part So = — 1 
which corresponds to non-interacting fermions. One writes 

S(p 2 , Pl ) = -1 + 2A(1 - e*™) du e -«(e*l+e-^-2Ae-i-* 2 )_ (4 g } 

Jo 

and integrates both pi and p 2 from to 2tt along the real axis before inte- 
grating over u. Both this definition of the integration and the choice for the 
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amplitude A(l\, I2) ensure that the initial condition is indeed satisfied in the 
physical domain k% > k\, I2 > h- With the constraint p\ + P2 = origi- 
nating in translational invarinace one recovers in this way the expression for 
the two-point correlator derived by Droz et al. |H[ by direct solution of the 
equations of motion with generating function techniques. 

Note that for the pure Glauber case A = and hence S — Sq — — 1. In 
this case the T> p anticommute and can represented e.g. by a Jordan- Wigner 
transformation of the usual local spin-1/2 raising or lowering operators s k . 
We recover in a surprising way the free fermion nature of this process. How- 
ever, since Glauber dynamics are well-understood []7T| we do not further 



pursue this matrix representation. In the present context we are more in- 
terested in the observation that for large times only small pi,p2 contribute 
to the integral (|4]). Making a substitution of variables Pi — > Pi/Vt and 
expanding for large t leads to 

s = -l + o(r 1/2 ) (44) 



for the late time behaviour of the scattering amplitude S (f43|). This proves 
our assertion for the two-point function: For large times the correction to the 
correlator due to coupling to the infinite-temperature Kawasaki heat bath is 
of subleading order 0(t -1 / 2 ) for any finite coupling £. 

Higher order correlation functions are treated analogously. The permu- 
tation of T> p matrices using the matrix relations (^), ( p3|) yields the Bethe 
ansatz equations and functional equations for the matrix elements T({pi}). 
For a r-point correlator this functional equation is solved by the r-particle 
Bethe wave function |4"0fl . Because of the integrability the scattering 
amplitudes factorize into products of two-body amplitudes. Thus for large 
times the leading part comes from the free- fermion amplitude S = — 1 and 
therefore the leading part of the correlator is independent of the Kawasaki 
coupling ( and given by the pure Glauber correlation function. The leading 
correction which results from this coupling is of order t -1 / 2 . This proves our 
assertion. 

5 Conclusions 

A 10-parameter class of stochastic reaction-diffusion systems can be mapped 
by a similarity transformation to a generalized Heisenberg quantum chain, 
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the spectrum of which can be obtained by the Bethe ansatz. It turns out 
that time-dependent r-point density correlation functions are given by the 
/ < r-magnon sectors. As an alternative to that approach we introduced 
a dynamic matrix ansatz for the general 12-parameter model. This ansatz 
reduces the calculation of all correlators to the calculation of certain matrix 
elements. These matrices satisfy an infinite-dimensional algebra which is 
determined by the bulk dynamics of the process. The boundary conditions 
determine which matrix elements one has to take. 

In the cases of the symmetric exclusion process and of the non-equilibrium 
spin relaxation model of Droz et al. the algebra satisfied by the matrices 
can be used to obtain the spectrum of the corresponding quantum Hamil- 
tonian. We constructed a two-dimensional time-dependent matrix represen- 
tation of the algebra for the symmetric exclusion process from which one 
obtains explicitly all r-point density correlators for a one-parameter class of 
initial states. The corresponding eigenvectors of the Heisenberg chain are the 
ground state with energy and a bound state with energy e = (a+P+ r y+5)/2. 
An alternative treatment of the dynamical algebra exploits directly the alge- 
braic relations which result in functional equations for the dynamical part of 
the correlator. We have solved these equations for the spin relaxation model 
in terms of Bethe wave functions. We proved that independently of the cou- 
pling to the infinite-temperature heat bath all r-point equal-time correlation 
functions decay to leading order in time like the zero-temperature Glauber 
correlators. 

^From a mathematical point of view the dynamical matrix algebra and 
its representation theory is not yet well-understood. The stationary version 
of the extended algebra was considered by Hinrichsen et al. |77j who con- 
structed a four-dimensional representation for a coagulation/decoagulation 



model. Krebs and Sandow job] could prove that the stationary algebra ex- 
tended with the auxiliary matrices forms an equivalent formulation of the 
stationary master equation. This guarantees the existence of a representa- 
tion for the general case. To date there is no equivalent theorem for the dy- 
namical algebra and no representation theory. A second important question 
concerns the relationship between the integrability of quantum chains and 
the dynamic matrix ansatz which emerged here and in other work |]2], 59, (Tgl - 
The result of Krebs and Sandow for the stationary algebra shows that there 
is no general relationship between the possibility of an algebraic description 
and the integrability of a system. Yet for integrable models the algebra is 
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powerful enough not only to recover the known Bethe ansatz equations but 
also to obtain results which cannot be obtained using standard Bethe ansatz 
techniques. This suggests that matrix algebras describing integrable models 
have special, as yet undiscovered properties. 
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